Magnetic trapping of neutrons 
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I Accurate measurement of the lifetime of the neutron (which is unstable to beta 

' decay) is important for understanding the weak nuclear forceB and the creation of 

T— I . matter during the Big Bangjl. Previous measurements of the neutron lifetime have 

O ; 

O I mainly been limited by certain systematic errors; however, these could in principle 

o . 

^ ' be avoided by performing measurements on neutrons stored in a magnetic traj 

' I ■ Neutral and charged particle traps are widely used tool for studying both composite 

Q ; 

^ . and elementary particles, because they allow long interaction times and isolation 

^ I from perturbing environment^ . Here we report the magnetic trapping of neutrons. 

^ ' The trapping region is filled with superfluid ^He, which is used to load neutrons 

H : 

- - ' into the trap and as a scintillator to detect their decay. Neutrons have a lifetime 

in the trap of 750^200 seconds, mainly limited by their beta decay rather than 

trap losses. Our experiment verifies theoretical predictions regarding the loading 

process and magnetic trapping of neutrons. Further refinement of this method 

should lead to improved precision in the neutron lifetime measurement. 

A more precisely known value of the neutron lifetime will benefit both weak interaction 

and Big Bang Nucleosynthesis (BBN) theory. Combined with the measured neutron beta-decay 
"''Present address: Institut Laue-Langevin, B.P. 156-6 rue Jules Horowitz, 38042, Grenoble (Cedex 9), France. 
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asymmetry, the neutron hfetime yields the Vud element of the Cabibbo-Kobayashi-Maskawa 
(CKM) quark mixing matrix. Non-unitarity of the CKM matrix would indicate physics beyond 
the Standard ModelB. In addition, the theoretical uncertainty in the predicted ^He abundance 
in the BBN model is dominated by the uncertainty in the neutron lifetimJi. 

Many experiments, using different methods, have measured the neutron beta-decay 
lifetimefi-i. The most precise measurements use ultracold neutrons (UCN, neutrons with ener- 
gies E/kh ^ 1 mK, where fcf, is Boltzmann's constant )0 confined in a 'bottle' made of materials 
with high UCN reflectivity. The largest uncertainties are due to UCN interactions with the 
material walls and fluctuations in the initial filling densitiesBi. Extrapolation methods are nec- 
essary to extract the neutron beta-decay lifetime. In another experiment, neutrons with specific 
trajectories and in a band of velocities (typically between 10-14 m s~^) were kept in a magnetic 
storage ring similar to those used in high-energy particle physicsi. Since storage rings only 
confine neutrons in two dimensions, coupling between the large longitudinal momentum and 
small radial momentum by way of betatron oscillations allows neutrons to escape. In contrast, 
neutrons in magnetic traps are not susceptible to wall losses or betatron oscillations. 

Static magnetic traps are formed by creating a magnetic field minimum in free space. 
The confining potential depth (D) of such a trap is determined by the magnetic moment of the 
trapped species (/i) and the difference (AB) between the magnitude of the field at the edge of 
the trap and at the minimum, D = /lAB. A particle in a low- field-seeking state (one with its 
magnetic moment anti-parallel with the local magnetic field vector) is pushed towards the trap 
minimum. Low-field-seeking particles with total energy less than D are energetically forbidden 
from leaving the trapping region. For atoms and molecules with a magnetic moment of one 
Bohr magneton it is possible to produce trap depths of ~ 1 K. The trap depth for a neutron in 
the same trap would be only ~ 1 mK, because of its much smaller magnetic moment. Despite 
this difficulty, magnetic trapping of the neutron was proposed as early as 1961 by VladimirskiS. 
His proposed technique was later used to confine neutrons using a combination of gravity and 
magnetsEi. A separate effort to trap neutrons using a similar loading method to our work (but 
different detection scheme) was unsuccessful because of the high temperature of the helium 
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during the loading phasdll. 

Crucial to the utility of traps are the techniques used to load them. In order to catch 
a particle in a static conservative trap, its energy must be lowered while it is in the potential 
well. Atoms and molecules have been cooled and loaded into magnetic traps by scattering 
with either cryogenic surfacesQIli, cold gaseJll or photons from a laser beam (laser cooling)0. 
Neutrons, however, cannot be loaded by such methods because they cannot be excited optically 
and interact too weakly with atoms to be effectively cooled by a gas. Direct thermalization 
with a cold solid or liquid is generally precluded by the high probability for neutron absorption 
in the vast majority of materials. 

Our trapping of neutrons relies on a loading technique that employs the "superthermal 
process"0. A neutron with kinetic energy near 11 K (where the free neutron and superfluid 
helium dispersion curves cross) that passes through the helium-filled trapping region can lose 
nearly all of its energy through the creation of a single phonon. Neutrons that scatter to energies 
less than the trap depth (^ 1 mK) and in the appropriate spin state are trapped. Upscattering 
of UCN in liquid helium by single phonon absorption is highly suppressed due to the low density 
of 11 K phonons in the superfluid (proportional to e"^^^/"^). The rate of higher-order (multiple- 
phonon) upscattering processes is proportional to T''(ref. |l^). The helium is maintained below 
a temperature of 250 mK to minimize upscattering of the trapped neutrons (total calculated 
upscattering rate per neutron < 10~^ 

Snother possible loss mechanism for trapped neutrons is nuclear absorption due to '^He 
contamination of the superfluid ^He bath. Unlike ^He which does not absorb thermal neutrons 
at all, '^He readily absorbs neutrons with a thermal cross section, a = 5330 x 10~^^ cm^. To 
minimize these losses, our superfluid-helium bath was isotopically purified (^He/'^He < 5x 10"^^; 
ref. ^), giving a calculated absorption loss rate per neutron of 5 x 10^^ s^^. Thus, the trapped 
neutrons travel unimpeded through the helium and reside in the trap for a time limited primarily 
by the beta-decay lifetime of the neutron. 

In our experiment, the superfluid helium is contained in a tube that is axially centered 
within a superconducting magnetic trap (see Fig. IT]). The incident cold neutron beam is colli- 
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mated by a neutron-absorbing ring, passes through the trapping region, and is absorbed by the 
beam stop. As the beam traverses the trapping region, about 1% of the 11 K neutrons scatter 
in the hehum. Some of these neutrons are trapped and the remainder are absorbed by shielding 
materials that surround the hehum. 

When a trapped neutron decays (into an electron, proton and anti- neutrino) , the result- 
ing high-energy electron (up to 782 keV) travels through the helium leaving a trail of ionized 
helium atoms. These ions quickly combine with neutral helium atoms to form excimer diatomic 
molecules, about half in singlet states and half in triplet states. About 15 singlet molecules 
are created per kiloelectronvolt of electron energy and the trail is 0-20 mm long0. The sin- 
glet molecules decay within 10 ns, emitting a pulse of light in the extreme ultraviolet (EUV), 
A ~ 70 — 90 nm. This pulse of scintillation light is the signal of a neutron decaying in the trap. 
Radiative decay of the triplet states has a much longer lifetime, which we have measured to be 
13 =b 2 s (ref. and does not contribute to our signal. 

The EUV scintillation light is wavelength-shifted to the visible and transported to the 
outside of the cryostat where it is detected. As shown in Fig. |l|, the superfluid-helium-filled 
trapping region is surrounded by an acrylic tube. The inside surface of this tube is coated with 
a thin layer of polystyrene doped with the organic fluor tetraphenyl butadiene (TPB; ref. ^3|) . 
The TPB converts the EUV scintillation light into blue light, a portion of which is internally 
reflected down the length of the tubeB. The tube is optically coupled to a solid light guide 
which transports the light to the end of the 250 mK region. The light then passes through a 
window at 4 K and into a second light guide which exits the dewar. The light is split into two 
guides that are each coupled to a photomultiplier tube (PMT). Background from uncorrelated 
photons is suppressed by requiring coincident detection of at least two photoelectrons in each 
PMT. This is required because of significant time-dependent luminescence (initial counting rate 
without coincidence detection of ~ 10^ induced by neutron absorption in the shielding 

materials. In a series of calibration runs (using a ^^^Sn, 360-keV beta line source placed in 
the center of the trapping region) we have determined that a total of 14.5 it 1.5 photoelectrons 
per megaelectronvolt of beta energy are detected in a pulse 20 ns wide. We estimate that a 
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coincidence signal is recorded for 31 it 4% of neutron decays in the trapping region. 

In each experimental run, the cold neutron beam is allowed to pass through the trapping 
region for 1,350 s, after which the beam is blocked and pulses of light are counted for 3,600 s. 
While the beam is on, neutrons interact with the helium and the number of UCN in the trap 
increases. After the beam is turned off, the PMTs are turned on and a decreasing rate of 
scintillation events is observed. A background signal, consisting of both constant and time- 
varying components, obscures the trapped neutron signal. These backgrounds are subtracted 
by taking data in two different kinds of run, trapping (signal + background) and non-trapping 
(background only), and subtracting the two. In a trapping run the magnetic field is on during 
the initial loading phase, so that UCN are confined in the trap, while in a non-trapping run 
the magnetic field is initially off, so that no neutrons can be trapped. In a non-trapping run, 
the field is turned on during detection so that any residual luminescence (which depends on 
the magnetic field) will be properly subtracted. Equal numbers of each kind of run were taken, 
pooled and subtracted to give the background subtracted data. 

Two sets of background subtracted data were collected: set 1 with a trap depth of 
0.76 mK (figure |2|a) and set 2 with a lower trap depth of 0.50 mK (figure ^). (The lower 
trap depth was used due to problems with the magnet.) Most of the run-to-run variation in 
background rate is eliminated by excluding the first two pairs of runs in which the background 
rate is changing quickly due to activation of materials with lifetimes > 12 h. The remaining 
23 pairs in set 1 (from about five days of running) and 120 pairs in set 2 (from about three 
weeks of running) are pooled and modeled as: 

VFi = Qie^*/" + Ci, W2 = a2e-*/" + C72. (1) 

The subscripts 1 and 2 refer to sets 1 and 2, ai = eNi/r, Ni is the initial number of trapped 
neutrons and i = 1 or 2, e is the detection efficiency, and r is the lifetime of neutrons in the 
trap. The constant C, is present to account for the possible remaining effect of the changing 
background rate due to long-lived activation. However, in all of our fits the value of Cj is 
consistent with zero. The fit is performed simultaneously on the two data sets, minimizing the 
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total while varying five parameters: ai, ai-, Ci, C2 and r. The only parameter connecting the 
two data sets is r. The best fit values indicate A^i = 560 =b 160 and N2 = 240 it 65. Calculations 
using the known beam flux, trap geometry and the theory of the superthermal process predict 
A'^i = 480 ± 100 and = 255 it 50, in good agreement with the measured values. The best 
fit value for the lifetime, r = 750^200 

s is consistent with the presently accepted value of the 
neutron beta-decay lifetime of 886.7 ± 1.9 s (ref. |). All of the errors quoted correspond to a 
68% confidence interval. 

To verify that our signal is due to trapped neutrons, we doped the isotopically pure ^He 
with ^He at a concentration of 2 x 10"'' '^He/^He. This amount of ^He absorbs the trapped 
neutrons in less than 1 s without affecting anything else in the experiment (less than 1% of the 
11 K neutrons are absorbed by '^He). The ^He background subtracted data (see Fig. ^3 and d) 
is consistent with zero, and differs from the signal data by more than four standard deviations. 
This confirms that our signal is due to trapped neutrons. 

This work demonstrates the loading, trapping, and detection techniques necessary for 
performing a neutron lifetime measurement using magnetically trapped UCN. Another impor- 
tant result is the direct confirmation of the theoretical prediction of the UCN density in the 
trap. Our value for the number of trapped neutrons at a 0.76 mK trap depth corresponds to a 
density of 2 UCN/cm^, compared to the density of 1 UCN/cm'^ obtained in the UCN material 
bottle experimentsS'i with a comparable UCN cut-off energy and higher flux reactor. Our mea- 
sured density is consistant with previous measurements of the UCN production rate based on 
observation of upscattered UCN at higher temperaturesH and confirms that the low number of 
UCN detected in ref. |2^ was not due to an intrinsically low UCN production rate. 

Magnetic trapping should allow significant improvement in the measurement of the neu- 
tron lifetime Tn- There are several loss mechanisms specific to this method, including thermal 
upscattering and '^He absorption. Both of these can be suppressed to a high level as discussed 
above. Another possible loss mechanism is neutron depolarization. When the neutron passes 
through a region where the field changes direction quickly compared to the neutron's spin pre- 
cession frequency, the neutron can flip spin and be ejected from the trap. In our trap geometry. 
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the trap has no zero-field regions and this loss rate can be minimized by the application of a 
bias fielJ i. 

Using our method with a higher-flux source of cold neutrons, a 10^^ measurement of 
r„, should be possible!. In the near term, improvements such as increased trap size and depth 
(resluting in an improved ratio of signal-to-background events) should allow an improved ac- 
curacy measurement of t„ at the 5 x 10~^ level. Other potential applications of superthermal 
UCN sources include searches for time-reversal violationH and new methods of neutron scat- 
tering especially suited for the study of large (biological-scale) moleculesS. 
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Figure 1: Half section view of the neutron trapping apparatus. The trapping region is filled 
with isotopically pure ^He at a temperature ^ 250 mK. The helium is contained within a 
cupronickel tube. A beam of cold neutrons passes through a series of teflon windows0 and 
enters the helium from the left. It is collimated by a boron carbide ring, passes through the 
trapping region and is absorbed by a boron carbide beam stop. Approximately 1% of the 11 K 
neutrons scatter in the superfluid helium. Those neutrons (yellow) in the low-field-seeking spin 
state and with energy below the trap depth are magnetically confined. The rest of the scattered 
neutrons are absorbed by neutron shielding material (boron nitride) surrounding the trapping 
region. The magnetic trapping field is created by an assembly of superconducting magnets. 
Radial confinement is provided by a quadrupole constructed from four racetrack-shaped coils 
and axial confinement is provided by two sets of solenoids. Electrons from neutron beta-decay 
cause EUV scintillations in the superfluid helium which are wavelength shifted to the visible 
by a thin film of TPB-dopedEl polystyrene coated on the inside of an acrylic tube surrounding 
the trapping region. This tube is optically coupled to an acrylic lightguide which transports 
the blue light to the end of the 250 mK region (to the right). 
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Figure 2: Counting rate as a function of time af- 
ter the neutron beam is turned off (pooled back- 
ground subtracted data). The two sets of trap- 
ping data are fit for five parameters: cti, a2, Ci, 
C2 and r. The constant offset terms, Ci and C2 
arc found to be -0.004±0.015 and -0.006zb0.006 
respectively, both consistent with zero. The best 
fit values for the initial counting rates for trapping 
data sets 1 and 2 are ai = 0.196 it 0.05 and 
02 = 0.084 ±0.02 s-^ The detection efficiency of 
e = 31 lb 4% combined with the known neutron 
lifetime, gives Ni = 560 ±160 and N2 = 240 ±65. 
The best fit value for the Hfetime is r = 750^200 s- 
The results of this fit are shown on plots a and c. 
The •^He data is fit to four parameters, ai, 02, Ci, 
C2, with T fixed at 750 s. Because no fit parameter 
applies to both data sets, this is equivalent to two 
separate two-parameter fits, one to each data set. 
The amplitudes, af ^ = 0.040 ± 0.045 and ^ = 
0.005 ± 0.016, are both consistent with zero. The 
constant offset terms are Ci = 0.011 ± 0.011 and 
C2 = -0.006 ± 0.004. The ^He data and fit are 
shown on plots b and d. (a) Trapping data set 1 
(23 pairs of runs at a trap depth of 0.76 mK). (b) 
^He data set 1 (20 pairs at 0.76 mK). (c) Trap- 
ping data set 2 (120 pairs at 0.50 mK). (d) ^He 
data set 2 (125 pairs at 0.50 mK). 
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